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Abstract — A Generalised Hadamard Transform for multi- 
phase or multilevel signals is introduced, which includes the 
Fourier, Generalised, Discrete Fourier, Walsh-Hadamard and 
Reverse Jacket Transforms. The jacket construction is formalised 
and shown to admit a tensor product decomposition. Primary 
matrices under this decomposition are identified. New examples 
of primary jacket matrices of orders 8 and 12 are presented. 

I. Introduction 

There are many discrete signal transforms whose transform 
matrices have entries on the unit circle. For instance, the 
Discrete Fourier Transform (DFT) for signals of length n 
and the Walsh-Hadamard Transform (WHT) are both special 
cases of the Fourier Transform (FT) found by interpreting 
the Cooley-Tukey Fast Fourier Transform in terms of abelian 
group characters [6]. The family of discrete Generalised Trans- 
forms {(GT) r ,0 < r < m - 1} for signals of length 2 m [3, 
10.2] includes the WHT as case r = 0, the complex BIFORE 
transform (CBT) as case r = 1 and the 2 m -point DFT as case 
r = m — 1. Both the WHT and DFT are suboptimal discrete 
orthogonal transforms, but each has wide application. 

These transform matrices all fall into the class of Butson 
Hadamard matrices [2]. 

For signals where entries outside the complex n th roots 
of unity are needed, Lee has introduced multi-phase or mul- 
tilevel generalisations of the WHT and of the even-length 
DFT under the name Reverse Jacket Transforms (RJT), see 
[7]. These are so-called because the border ("jacket") and 
the centre of the unitary matrix representing the transform 
change independently under inversion. The former admit a 
recursive factorisation into tensor products so represent a fast 
transform similar to that of the WHT. (Note that the formula 
[7, Definition 5] is not a reverse jacket transform generalising 
both the WHT and even-length DFT, as claimed, because it is 
not unitary.) 

The class of Butson Hadamard matrices can be generalised 
to include all these transform matrices. 

II. A Generalised Hadamard Transform 

In this most general situation, we work in a ring R with 
unity 1. This includes R, C and Galois Field alphabets 
GF(p a ), though if we need to distinguish signal values x and 
—x, the ring must have characteristic ^ 2. 

Definition 2.1: Suppose R is a ring with unity 1, group of 
units R* and that char R does not divide v. A square matrix 



M of order v > 2, with entries from a subgroup A < R* is 
a Generalised Butson Hadamard ( GBH) matrix, if 

MM* = M*M = vl v , 

where M* is the transpose of the matrix of inverse elements of 
M: m*j = (mj,) -1 . It is denoted GBH{N, v), or GBH(w, v) 
if A is finite of order w. 

A GBH matrix is always equivalent to a normalised GBH 
matrix, which has first row and column consisting of all Is. By 
taking the inner product of any non-initial row of a normalised 
GBH matrix M with the all- Is first column of M*, we see 
that the sum of the entries in any row of M, apart from the 
first, must equal 0, and similarly for rows of M* (columns of 
the matrix of inverses M^ -1 ' = [m^- ]). The tensor product 
of two GBH matrices over the same group A is a GBH matrix 
over A. 

Definition 2.2: Let x be a signal of length n from R*, 
where n G R*, let A < R* and let B be a GBH(N, n). 
A Generalised Hadamard Transform (GHT) of x is 



x = B x 



(1) 



and an Inverse Generalised Hadamard Transform (IGHT) of 
x is 

x = n" 1 S*x. (2) 
The next section describes a construction for GBH matrices 
of even order and additional internal structure, containing the 
WHT, even-length DFT and reverse jacket transform matrices. 

III. The jacket matrix construction 

Throughout this section, let G be an indexing set of even 
order 2n (sometimes G is a group such as Z2„ or ZJJ but G 
may be non-abelian or the group structure may be irrelevant). 

Definition 3.1: Let R be a ring with unity 1. A normalised 
GBH(A,2n) matrix K indexed by G = {l,...,2n} with 
entries from A < R* is a jacket matrix if it is of the form 



K = 



1 ±1 



1 1 

* ±1 

* ±1 
±1 ±1 



(3) 



where the central entries * are from A. The jacket width 
of K is m > 2 if K is permutation equivalent to a jacket 
matrix K in which rows 2, . . . , m, 2n — m + 1, . . . , 2n and 



columns 2, . . . , m, 2n — m + 1, . . . , 2n all consist of ±1 and 
to is maximal for this property. Otherwise the jacket width is 

to = 1. 

All non-initial ±1 rows and columns of a jacket matrix if 
necessarily sum to in R. If if is a width to jacket matrix 
it follows that K* is itself a width m jacket matrix. If also 
2n e R*, then K has an inverse K^ 1 = (2n)^ 1 K* over R. 
If (2n)i £ R* then £ = (2n)-^K is unitary (^£T* = i 2 „) 
and we will (by a slight abuse of terminology) also say K is 
unitary. 

Example 3.2: The matrix St of the WHT of length 2* 



Si = 



5* = ®*5i, t>2, 



(4) 



is an extremal case (jacket width 2* 1 ) of such jacket matrices 
("all jacket"). 

Example 3.3: The matrix C t of the CBT of length 2* 



Ci = 



C — 



1 -', 
1 \ 

Si Si 
Ci -Ci 

Cm — 1 

Ci ® <S m 



(5) 



-Ci 



,to > 3. 



is permutation equivalent to a jacket matrix, for t > 2. 
Proof. For t > 2, Ct is normalised and the first two rows of 
C t consist of 2 t_1 copies of Si. By induction the (2* _1 + l) s * 
column of C t is [1 -1] T , where 1 has length 2* _1 . Rotating 
the second row to the bottom and the (2* _1 + l) st column to 
the right of C t produces a jacket matrix. □ 

Example 3.4: The matrix Tin = [ ^ k ]o<j,fe<2n-i> where 
u> = e -7 ™/™, n > 1, of the DFT of length 2n is permutation 
equivalent to a jacket matrix. 

Proof. [7, Theorem 1, Definition 1] Represent the indices in 
mixed radix notation j = j\n + jo = (ji,jo), with index set 
G = Z 2 x Z n . The permutation (ji, j ) >-> (ji, (1 - jt)j + 
(n—l—jo)ji) leaves the first n indices unchanged and reverses 
the order of the last n indices. Under this permutation on rows 
and columns, Ti n is equivalent to )C n (ui) = 

J w {(l-Ji)j'o+(n— l-joMi+jin}{(l-fci)fco + (n-l-fco)fci+fcin}j 

(6) 

the matrix of Lee's complex RJT □ 

If K, K' are jacket matrices indexed by G, G' of orders 
2n, 2n' respectively, with entries from R*, then the tensor 
product K ® K ' is a jacket matrix indexed by G x G' of 
order Ann', with entries from R*, since the border condition 
is easily seen to be satisfied. In fact, a tensor product of jacket 
matrices is a jacket matrix of width > 2. 

Lemma 3.5: If Ki is a width TOj jacket matrix with entries 
from R*, for i = 1,2, then ifi <g> if 2 is a jacket matrix of 
width at least 2toito 2 . 

Proof. Let ifi have order 2n and if 2 have order 2n'^ermute 
Ki to ifi, so K\ ® Ki is permutation equivalent to K\ <g> K2. 
Let i G {1, . . . , mi, 2n — mi + 1, . . . , 2n} be an index of 
an all-(±l)s row in K\. The corresponding i th block row in 



K\®Ki consists of 2n copies of ±lif 2 , so each row indexed 
{2, . . . , m 2 , 2n' — m 2 + 1, . . . 2n'} of each copy consists of 
all-(±l)j^contributing 2m 2 — 1 all— (±1) rows to the i th block 
row of K\ ® if 2 . If i = 1 the top row is all Is and if i > 1 
the top row consists of n Is and n —Is. Those in the top mi 
block rows of K\ ® if 2 may be permuted to occupy the top 
2toito 2 rows and those in the bottom mi block rows to the 
bottom 2mim 2 rows, and similarly for columns. □ 

IV. Tensor decomposition 

If a jacket matrix may be decomposed as a tensor product 
of two smaller jacket matrices, the decomposition may be 
repeated until no further tensor product decomposition is 
possible. 

Definition 4.1: A jacket matrix of length 2n is a primary 
jacket matrix K n if it is minimal with respect to tensor 
product, that is, there are no jacket matrices K, K' such that 
K n is permutation equivalent to K <g> K' . 

Examples of primary jacket matrices for n = 1, . . . , 4 are 

K x =Si, 
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r ±1 e R* (7) 



(8) 



where a is a primitive 6 th root of unity in an integral domain 
R, and 
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(9) 



The matrix K\ is the unique 2x2 jacket matrix, and 
if 2 (l) = K\ ®K\ — 5 2 , so is not primary. The matrix if 2 (r) 
for r ^ ±1 e M* is a "centre-weighted Hadamard transform" 
(CWHT) matrix [7] and for r = i e C*, is K 2 (i). The matrix 
K 3 (a) with a = e n l z is /C 3 (e™/ 3 ) and with a the fourth 
power of a primitive root in GF(25) is an "extended" complex 
RJT matrix [7, Example 4]. The matrix K±{i) is described for 
the first time here: it is equivalent to the back-circulant matrix 
derived from a quadriphase perfect sequence [1, Example 2] 
of length 8; such sequences are rare. 

Corollary 4.2: A jacket matrix is permutation equivalent 
to a tensor product of one or more primary jacket matrices. 



Conversely, any tensor product of primary jacket matrices is 
a jacket matrix. 

V. Construction of primary jacket matrices 



By Lemma 13751 any jacket matrix of width 1 is primary. The 
tensor product of two jacket matrices is a jacket matrix, but in 
fact, to construct a jacket matrix it is enough that one factor 
is a jacket matrix and the other a normalised GBH matrix. 

Theorem 5.1: Let B be a normalised GBH matrix and K 
a jacket matrix, both with entries in R*. Then B ® K is 
permutation equivalent to a jacket matrix (B <g) Ky. 
Proof. Let B have order m and K have order 2n. The first 
2n rows of B <g> K consist of m blocks K, K, . . . , K, so the 
I s * row is all l's and the 2n th row is all ±ls, and similarly 
for columns. Cyclically permute row 2n to row 2mn and row 
2n+i to row 2n+i— 1, i = 1, . . . , 2n(m— 1), and similarly for 
columns. This shifts row 2n to the bottom of the matrix and 
column 2n to the right of the matrix, leaving the order of the 
other rows and columns otherwise unchanged. This permuted 
matrix (B <g) K)^ is of the form @. □ 

This result explains the generation of some primary jacket 
matrices and is fundamental to the construction of Generalised 
Hadamard Transforms. 

Example 5.2: Let (3 f 1 G R* satisfy f3 2 + (3 + 1 = 0, so 
/3 3 = 1. Let 
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, 3). Then 
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This jacket matrix relates to the DFT matrix of (|8} as fol- 
lows. A second permutation (2543) cycling central rows and 
columns gives the jacket matrix 
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When a is a primitive 6 th root of unity in R* with a 2 — p 
and a 5 = 7, then this matrix equals £3(7). 

The following matrix is a new 12x12 primary jacket matrix, 
since it has width 1. 

Example 5.3: Let r^ile R*. 

Let K 6 (p,r) = (B 3 ®K 2 (r))i = 



a 

-7-0 

r/3 
-0 

2 
r0 2 
r 2 

-0 2 
- 1 



0" 
r0 2 
r0 2 
-0 2 


-T0 
T0 
-0 



0' 
T0 2 
-T0 2 
~0 2 



Those GHT matrices which are jacket matrices have addi- 
tional structure, by virtue of their tensor product decomposi- 
tion into primary jacket matrices and their jacket form, which 
may particularly suit them to specific applications. 

Consider the set of jacket matrices with entries in in an 
integral domain R 

{K = (p'Kt) ® K 2 (rY ® IC n {a) s ; £>0,e,Se {0, 1}}, 

(10) 

where r ^ il e Ji*, a is a primitive 2n th root of unity 
and where by M° we mean the lxl identity matrix. When 
£ > 1, e = 0, 5 = 0, this is the WHT When £ = 0,e = 0,<5 = l 
and R = C this is equivalent to the 2n-point DFT. When 
e = 1, 5 = and R = R, this is the CWHT When e = 0, n = 
2,a = i € C*,5 = 1, or when e = 1, r = i G C*, 5 = 0, 
this is the complex RJT, and when e = 0, S = 1, this is the 
extended complex RJT. 

This uniform classification may make it possible to recog- 
nise common fast algorithms. 

Finally, the examples and constructions of primary jacket 
matrices above are all instances of a class of matrices called 
cocyclic, investigated by the author and colleagues over the 
past decade [4], [5]. 

VI. Conclusion 

To summarise: a generalisation of Butson's Hadamard ma- 
trices determines a Generalised Hadamard Transform (GHT). 
The GHT includes the Fourier and Generalised Transform 
families (in particular the WHT and DFT) and the centre- 
weighted Walsh-Hadamard, Complex Reverse Jacket and ex- 
tended Complex Reverse Jacket Transforms. 

In the jacket case, GHT matrices can be permuted into 
tensor products of primary jacket matrices. New primary 
jacket matrices may be constructed as tensor products of a 
Generalised Butson Hadamard matrix which is not a primary 
jacket matrix, and a primary jacket matrix. New examples in 
orders 8 and 12 have been given. 

Application of the GHT to image processing, error-control 
coding and decoding and sequence design are obvious direc- 
tions for future research. 
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